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Abstract 

Recently, Jadbabaie, Lin, and Morse (2003) offered a mathematical 
analysis of the discrete time model of groups of mobile autonomous agents 
raised by Vicsek et al. in 1995. In their paper, Jadbabaie et al. showed 
that all agents shall move in the same heading, provided that these agents 
are periodically linked together. This paper sharpens this result by showing 
that coordination will be reached under a very weak condition that requires 
all agents are finally linked together. This condition is also strictly weaker 
than the one Jadbabaie et al. desired. 

Index Terms — Decentralized control, multi-agent coordination, switched 
systems. 



1 Introduction 

Coordination of groups of mobile autonomous agents (particles |9., or boids jS]) 
has attracted researchers in a surprisingly wide variety of disciplines ranging 
from physics 0011], to the biological sciences [HUE], to computer science and 
engineering El EH Hj • 

This paper is mainly concerned with one particular discrete time model of 
groups of mobile autonomous agents, viz., the one proposed by Vicsek et al. [H] 

*This work was partly supported by the National Foundation of Natural Science of China 
(60305005, 60321002, 60496321) and by a Hong Kong CERG (Grant No. CityU 1234/03E). 
I 1 Email: lisanjiang@tsinghua.edu.cn (S. Li), iswang@cityu.edu. hk (H. Wang) 
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in 1995. In this model, a group of autonomous agents is moving in the plane with 
all agents moving at the same speed but with different headings. Each agent's 
motion is updated using a local rule based on the average of its own heading and 
the heading of its "neighbors." This is known as the nearest-neighbor rule in [3]. 
Agent z's neighbors at time t are those agents that are either in or on a circle of 
pre-specified radius r centered at agent z's current position. Known as the Vicsek 
model, this can be viewed as a special version of a model proposed by Reynolds [5] 
for simulating animal aggregation for the computer animation industry. Although 
the Vicsek model is very simple, simulation results in [Oj show that, using the 
local update rule, all agents shall eventually move in the same direction despite 
the absence of centralized coordination, and that neighborhood of each agent will 
change. 

Recently, Jadbabaie, Lin, and Morse |3] offered a mathematical analysis of 
this model and provided a theoretical explanation for the observed behavior. 
They adopt a more conservative approach, which ignores how the neighbor-graphs 
depend on the agent positions in the plane. Note that under this assumption, 
the Vicsek model is a graphic example of a switched linear system. Their goal in 
that paper was to determine for a certain large class of switching signals and for 
any initial set of agent headings that the headings of all agents will converge into 
the same steady heading. 

Jadbabaie et al. jH] established sufficient conditions given in terms of neighbor- 
graphs for coordination of agents. One main result of jHj shows that all agents 
shall eventually move in the same heading if these graph are periodically jointly 
connected, i.e., the union of any T sequential graphs is connected for some fixed 
T. This is a nice result, but as Jadbabaie et al. put it (3J p990, below Theorem 
2], what one would prefer instead is to show that coordination would be reached 
eventually for every switching signal for which there is an infinite sequence of 
bounded, non-overlapping (but not necessarily contiguous) intervals across which 
the agents are linked together. 

This paper will show above condition desired by Jadbabaie et al. is indeed 
a sufficient condition for asymptotic convergence. This, however, follows from a 
more general observation: our main result in this paper shows that convergence 
will be attained if these neighbor-graphs are finally jointly connected, i.e., the 
union of all graphs started from any time is connected. 

The structure of this paper is as follows. In Section II we give a formal descrip- 
tion of the Vicsek model in terms of switching signals. Section III provides the 
major results, for both leaderless coordination and leader-following coordination. 
Conclusions and future work are given in the last section. 

2 The Vicsek model and the nearest-neighbor 
rule 

In this section, we review some basic definitions concerning the Vicsek model. 
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The system studied by Vicsek et al. [9 j consists of n autonomous agents, 
e.g., particles, robots, etc., labeled 1 through n. All agents move in the plane 
with the same speed but with different headings. The system operates at discrete 
time t — 0, 1, 2, • ■ • Let r > and v > be given numbers associated with the 
system. The dynamics of agent i is described by the sequence {xi(t),yi(t), 9i(t)}, 
where Xi(t),yi(t) G R are the coordinates of the agent in the plane, and 9i(t) is 
its heading taking value from [0, 2ir). At any time t = 0, 1, 2, • • ■ , each agent's 
heading is updated using a simple rule based on the average of its own heading 
plus the headings of its neighbors. For any two agents i,j, we say j is a neighbor 
of % at time t, written j G Ni(t), if d(i,j) < r, 1 where 

d(i,j) = yW) - xtf)) 2 + fait) - VM 2 - (!) 
Then, agent z's next heading is defined as 

W + 1 ) = TTW) () 

where rii(t) is the number of agent i's neighbors at time t. Agent i's next coordi- 
nates are defined as 

Xi (t+1) = x(t)+Vi(t)coa(9i(t)) (3) 
Vi(t + 1) = y(t) + Vi(t)an{9i(t)). (4) 

For any time t > 0, we define the neighbor- graph of the system described 
above as the simple undirected graph G(t) over V = {1,2, ••• ,n} where the 
vertex % corresponds to agent i and two vertexes, are connected by an edge 
in the graph if they are neighbors at time t, i.e., if j G J\fi(t). Since the neighbor 
relation can change over time, so can the graph that describes them. In the sequel, 
we write V for the collection of simple undirected graphs over V. A switching 
signal is a function a : N — > V that assigns to each time t a neighbor-graph that 
specifies the neighbor relation between agents. Clearly, for a Vicsek model, the 
function that assigns to each time t the neighbor-graph G(t) is a switching signal. 

Note that for the Vicsek model, the neighbor-graph is determined by the 
initial positions and headings of all agents as well as the pre-specified r > and 
v > 0. A complete description of the model would have to explain explicitly how 
a changes over time. As it is difficult to take this into account in a convergence 
analysis, Jadbabaie et al. adopt a more conservative approach, "which ignores 
how a depends on the agent positions in the plane and assumes instead that o 
might be any switching signal in some suitably defined set of interests." [H] 

We in this paper follow this basic assumption and formalize the Vicsek model 
as follows: 



1 In certain situation, choosing open neighborhood would give rise to more desirable results, 
see 0. 
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Definition 2.1 (Vicsek model). Given n agents, labeled 1,2, ••• , n, moving 
in the plane at discrete time t G N, let V be all simple undirected graphs over 
V = {1, 2, • • • , n}. A Vicsek model is a pair (V,a), where a : N — > "P is a 
switching signal. 

For each agent i, define i's a-neighborhood at time t, written J\fi(t), to be the 
set of agents that is connected to i by an edge in graph a(t). That is, agent j is 
a neighbor of agent z if and only if is an edge in graph u{t). 

Given an initial heading 9(0) = (6>j(0))™ =1 , agent z's heading 9i(t) evolves in 
discrete time according to Eq. El Namely, agent z's heading at time t + 1 is the 
average of the headings of agent % and its neighbors at time t. 

Remark 2.1. This definition of a multi- agent coordination model is very general 
and more flexible. Several dimensions of extension/completion could be incor- 
porated in this model: 1) We can choose either closed/open disk, a triangle-like 
zone or any subset of V as the neighborhood; 2) The velocity could also change 
in discrete time; 3) We could consider other state variables of agents besides 
their headings; and 4) The neighbor-graph could also be directed. This flexibility 
would be helpful in practical applications. 

The goal of this paper is to show for a large class of Vicsek models (or switching 
signals) and for any initial set of agent headings that the headings of all n agents 
will converge into the same heading. Compared with the results obtained in jH], 
ours are more general. 

3 A sufficient condition for multi-agent coordi- 
nation 

3.1 Notations and preliminaries 

Suppose (V, a) is a Vicsek model. Following Savkin 0, we define a graph cx(oo) 
over V — {1, 2, • • • , n} as follows: for any two nodes i, j, (i,j) is an edge in cr(oo) 
if and only if for any K > 0, there exists some k > K such that (i, j) is an edge in 
graph a(k). For convenience, given a collection of graphs {G x : x G X}, we write 
L±Lex ^ f° r the union of these graphs, i.e., any pair (i, j) is an edge in \+) x&x ^ 
if and only if it is an edge in some G x . Then it is easy to show that there exists 
some K > such that cr(oo) = \+) t>k cr(t) holds for all k > K. In this paper 
we will show that all agents shall eventually move in the same heading provided 
that cr(oo) is connected. This condition is more general than the one given in 
[Hj, where the authors require that the cr(t)'s are periodically jointly connected. 
In what follows, a switching signal a is called finally jointly connected if o~(oo) 
is connected. Clearly this is equivalent to saying that \+) t>k cr(t) is connected for 
any fcGN. 

For a sequence {/(&)} and a number u in R, we say u is an accumulation 
point of {f{k)} if there is a subsequence of {f(k)} that converges to u. We write 
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Accu({/(fc)}) for the set of accumulation points of {/(&)}. 

Given a Vicsek model (V, a) and an initial headings 0(0) = (0j(O))™ =1 , we now 
fix some notations concerning the model. 

For i — 1, 2, • • • , n, define 

9^ = {0i(t):*eN} (5) 
rrii = minAccuGj (6) 
Mi = maxAccuOj (7) 

n 

m = minmj (8) 

i=i 

M = max Mi. (9) 
i=i 

Note that 0j is a bounded set and, therefore, has minimum and maximum ele- 
ments. 



For any t 6 N, define 



0(t) = min 0i(i) (10) 



i=l 



0(t) = m&x6i(t). (11) 



i=l 



The following lemma shows 9_(t) < m < M < 6(t). 

Lemma 3.1. For anyt E N, we have 6 (t) < 9(t + l) < m < M < 9(t+l) < 9(t). 
Consequently, we have lim^ oo 0(t) = m and lim^oo 9(t) = M. 

Proof. For any non- negative t, note that by Vicsek's nearest-neighbor rule (Eq.|2J), 
we have 9(t) < 9 t (t+ 1) <J(t). In particular, we have 9(t) < 9(t + l) < 9(t + l) < 
9(t). Now since {9(t)} ({9(t)}, resp.) is a bounded ascending (descending, resp.) 
sequence, it has a limit. We now show its limit is m (M, resp.). Take {#_(£)} 
as an example. Since it is convergent, any subsequence of {9(t)} also converges 
to its limit. Suppose {/(&)} i s a sequence such that linn^oo 9i(f(k)) = m for 
some agent i. Note that 9(f(k)) < 9i(f(k)) for any k, we have lim^ oo 0(t) = 
lim^tx, 9(f(k)) < lim^oo 9i(f (k)) = m. On the other hand, since there exists 
some agent i such that {t : 9i(t) = 9(t)} is infinite, we have a sequence {g(k)} 
such that 9i(g(k)) = 9(g(k)). This shows that lim 4 ^ oo 0(t) = lim^oo 9(g (k)) = 
lim^^oo 9{(g(k)) > m since m is the minimum accumulation point. As a result, 
we have lim t ^ oo 0(t) = m. Similarly, we can show lim t __ 00 9(t) = M. So we have 
9(t) <m<M< 9(t) for any t. □ 

Note that as shown in the proof of the above lemma, we have a sequence, say 
{f(k)}, such that 9i(f(k)) = 9(f(k)) and lim^oo 9i(f(k)) = m for some agent i. 
Similarly, we have a sequence, say g(k), such that lim^oo 9 j(g(k)) = M for some 

3- 
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3.2 Leaderless coordination 

Theorem 3.1. Given a Vicsek model (V,cr), suppose o~(oo) is connected. Then 
for any 9(0) = (6i(0))2=i, we have 

lim 9i(t)=9 ss (i = l,2,... ,n) (12) 

t— >oo 

where 9 SS is a number depending only on 9(0) and a. 
To prove this theorem, we need several lemmas. 

Recall V = {1, 2, ■ ■ • , n}. For a graph G over V and any two disjoint subsets 
A, B of V, we say A and B are connected if there exist a G A, h G I? such that 
(a, 6) is an edge in G. If A happens to be a singleton {a}, we also say node a is 
connected to B. In this case, we say alternatively a has a neighbor in B. 

The following lemma suggests that, if the agents are divided into two parts 
such that the maximum heading of the first part is sufficiently smaller than the 
minimum of the second part, then, after updating the headings using Eq. the 
agents will also form two parts such that one part is still sufficiently smaller than 
the rest. 

For a < b in M and any natural number t, we write V t (a, b) = {i G V : a < 
9 l (t) <b}. 

Lemma 3.2. Given a < (3 < 7 and set 5 = [3 — a, e = 5/n n , suppose V t (a — 
e, a + e) and V t ((3 — e, 7 + e) are two nonempty disjoint subsets of V such that 
their union is V . We have 

(1) IfVt(a — e, a + e) and 14(/3 — e, 7+e) are disconnected at time t , thenVt+i(a — 
e, a + e) = V t (a - e, a + e) and V t+ i((3 - e, 7 + e) = V t (f3 - e, 7 + e). 

(2) If Vt(a — e, a + e) and V t ((3 — e, 7 + e) are connected at time t, then V t+ i(a — 
e, a + e) = V t (a — e, a + e) — {i G V : i has a neighbor in V t ((3 — e, 7 + e) at time t) 
and V — V t+ \(a — e, a + e) = V t+ \(a + 5/n — e, 7 + e). 

Proof. H V t (a — e,a + e) and Vj(/3 — e, 7 + e) are disconnected at time t, then for 
any i G V t (a — e, a + e), its neighbors are all in V t (a — e, a + e). By Eq. |2l we 
have a — e < 9, t (t + 1) < a + e. Similarly, for any j G V t ((3 — e, 7 + e), we have 
/3-e< 0j(* + l) < 7 + e- 

On the other hand, suppose V t (a — e,a + e) and Vt(/5 — e, 7 + e) are connected 
at time t. For i G V 4 (a — e, a + e), if all its neighbors are in V t (a — e, a + e), then 
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i G Vt+i(a — e, a + e); if i has a neighbor, say jo, in Vt(P — e > 7 + e ) 5 then we have 

0i(*) + £ i6 A/;(*)0i(t) 



> 



1 +Tli(t) 

^(t) + 9 jo (t) + ^j^sfiit^jjijoOjit) 
1 + rait) 

(a-e) + (p-e) + (m(t) - 1) x (a - e) 

1 + Tli(t) 

/3 — a + (1 + x (a — e) 



fa - e + 



l + n^t) 
/3 — a 

l + Tli(t) 

> a — e + S/n. 

Note that ^(t+l) < 7+e holds for any i & V. This shows that, if z G V t (a—e, a+e) 
has a neighbor in V*(/3 — e, 7 + e), then i G I4+i(a + #/ n ~ e ;7 + e )- Similarly, 
for any j G Vt(/3 — e, 7 + e), we can show 9j(t + 1) > a + — e. In summary, 
we have i G Vt+i(a — e, a + e) if and only if i G V t (a — e, a + e) and it has 
a neighbor in Vt(/3 — e, 7 + e) at time t. As for any other agent j, we have 
j G (a + 5/n — e, 7 + e) . □ 

Lemma 3.3. Suppose a (00) is connected and {/(&)} is a sequence. Then we 
have a subsequence {g{k)} of {/(&)} swc/i iZzcrf a// {#j(g(/c))} are convergent for 
ieV. 

Proof. This follows from the compactness of [0, 27r] and that G [0, 2ir) for 

any z, t. □ 

Lemma 3.4. Suppose {g{k)} is a sequence such that {9i(g(k))} converges to U for 
i — 1, 2, ■ • • , n. Then m = min™ =1 Zj and M = maxf =1 Zj, where m = lim^ 00 ^(t), 
M = hin^oo^t) and 9(t) = mi< =1 9i{t), 9{t) = max™ =1 ^(t). 

Proof. Take m = min™ =1 Zj as an example. Note that there exists some i such 
that {k : 9i(g(k)) = 9(g(k))} is infinite. We have a subsequence {/&(/>;)} of {g(k)} 
such that 9i{h{k)) = 9(h(k)) and Zj = lim^oo 9i(h(k)) = linn^oo 9 (h(k)) = m. 
That M — max™ =1 Zj is similar. D 

Proof of Theorem \S.l\ To show that these autonomous agents shall eventually 
move into the same heading, we need only to show m = M. We prove this by 
reduction to absurdity. 

Suppose m < M and {g{k)} is a sequence such that {9i(g(k))} converges to Zj 
for i — 1, 2, • • • , n. Recall by Lemma f3. 41 that m = min" =1 Zj and M = max™ =1 Zj. 
Set Z = min{Zj : Zj > m} and take 5 = I — m, e = 5/n n . Then there exists K > 
such that 9i(g(k)) G (Zj— e, Zj+e) for k > K and z = 1, 2, • ■ • , n. Moreover, we have 
V g (k)( m — e, m + e) = {i <E V : li = m} and V g ^)(l — e, M + e) = {i G V : Zj > I}. 
Clearly V g (k){m— e, m+e) and V^(fc)(Z— e, M+e) satisfy the condition of Lemma 13. 21 
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Now since i+J t>p cr(t) is connected for any p G N, we have some k > K such 
that V^(fc)(m — e, m + e) and V^(fc)(£ — e, M + e) are connected at time g(k)+w for 
some < w < g{k + 1) — g(k). Fix one such k and suppose g(k + 1) — #(&) = W. 
For each tu = 0, 1, • • ■ , W, we define A w = V g ( k ) +W {m — e, m + e) and B w = V — A w . 

Note that because Aw = V g ^+i)( m — e,m + e) = V g u.){m — e, m + e) = A , it 
is also true that m — e < 9(g(k)) < d(g(k + 1)) < m + e. Recall that because 6_{t) 
is an ascending chain (see Lemma I3.1J1 . we also have 9(g(k) +w) G (m — e, m + e) 
for any w — 1, 2, • • • , W — 1. 

Set C = {«)G [0, W 7 ) : v4 w and are connected at time + w}. Clearly 
C is not empty since there exists some w such that A Q is connected to B at time 
g{k) + w. Suppose C = {wi,W2, • ■ ■ , w g } and < w\ < w<i < ■ ■ ■ < w q < W. We 
claim 

(13) 

B Ws+1 = V g(k)+Ws+1 (m + — — e, M + e) (s — 1,2, - • ■ ,q) (14) 

As the induction basis, note that A = V g ^(m — e, m + e) = {i : U = m} and 
B = {i : h > 1} = V g(k) (l — e, M + e) = 7 p(fc )(m + 5 - e, M + e). 

Note that u>i is the first index 10 such that A w is connected to B w ; by 
Lemma 13 .2| we have A = A w and B w = V g ( k ) +W {m + 5 — e, M + e) for any 
w; < W\. Moreover, since A m = A is connected to B Wl = B at time g{k) + w%, 
by Lemma f3. 21 we have 

A Wl+ i = A m -{i6 A W1 : i has a neighbor in B Wl at time g(k) + wi} C A (15) 

= ^ - A Wl+1 = V g{k)+Wl+ i(m + 5/n - e,M + e) (16) 

Recall that A Wl+ i ^ since 6_(g(k) + w% + 1) G (m — e, m + e). 
Suppose for s < g we have 

(17) 

A^+i = V g{k)+Wj+1 {m + — - e, M + e) {j = 1, 2, • • • , s) (18) 
Note that s < n — 1 must hold since A contains at most n — 1 agents and 

A 2 Aoi+i 2 • ■ ■ 2 A«,+i ^ 0- 

We now show A Ws+x 2 A w a+1 +i and5 tos+1+1 = l / g (fc) +tos+1+ i(m+^ T -e, M+e). 

Note that w s+ i is the first index w > w s such that A w is connected to 5^. By 
Lemma El we have A w = A Wa+1 and B w = V g ( k)+W (m + ^ - e, M + e) = S Ws+1 
for any w G (w s , w s+ \]. Moreover, since A Ws+1 = A Wa+1 is connected to B Ws+1 = 
B Ws+ i at time g(k) + w s+ \, by Lemma I3~2l we have 

A i — A — 

{i G A Ws+l : i has a neighbor in B Wa+1 at time g(k) + w s+ i} C A Ws+ i (19) 

5 Ws+1+ i = V - A Ws+1+ i = V g(k)+Ws+1+1 (m + — r - e, M + e) (20) 
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In summary, we have obtained that A = A m D A Wq+ \. 

Note that if w q < W — 1, then A w and B w are disconnected for any w G 
(w q , W). By Lemma l3~2l again. we know A w = Aw for w G (w q , W]. In particular, 
we have A Wq+ i = Aw- On the other hand, if w q = W — 1, we also have A Wq+ i = 
Aw- 

This suggests that if m < M, then A ^ A W - This is a contradiction. So our 
assumption that m < M cannot hold. This ends the proof of this theorem. □ 

Remark 3.1. Note that if a : N — » V is a switching signal for which there exists 
an infinite sequence of bounded, non-overlapping (but not necessarily contiguous) 
intervals across which the n agents are linked together, then cr(oo) is connected. 
By the above theorem, we know all agents would eventually move in the same 
heading for this a. Consequently, this theorem shows the desired condition given 
in (3J p990, below Theorem 2] is a sufficient condition for asymptotic convergence. 

The hypothesis of Theorem 13-H however, is still not necessary. For example, 
if some a(t) is the complete graph over V, then a coordination could be achieved 
at time t + 1. But if it is not connected, a(oo) will have 1 < p < n connected 
components, say Gi, G2, • • • , G p . Similar to the argument given above for Theo- 
rem we can show for any h = 1, 2, • • • ,p, there exists a heading 6h such that 
lim^oo 6i(t) = 9 h for any i e G h . 



3.3 Leader-following coordination 

In |3] , Jadbabaie et al. also consider a modified version of Vicsek's discrete-time 
system, which consists of the same group of n agents as before except that one 
leader agent, labeled 0, is added. Agent moves at the same constant speed as its 
n followers but with a fixed heading 9q. Agent i then updates its heading using 
the average of its own heading plus the headings of its neighbors. Note that this 
time the leader may be in its neighborhood. 

Our abstract Vicsek model with a leader now can be formulated as follows: 

Definition 3.1. Suppose V + — {0, 1, • • • , n) and V + is the collection of simple 
undirected graphs over V + . A leader-following Vicsek model is just a pair (V + , a), 
where a : N — ► V + is a switching signal. 

For each agent % > 0, define i's cx-neighborhood at time t, written A/i(t), to be 
the set of agents that are connected to i by an edge in the graph a(t). That is, 
agent j is a neighbor of agent i if and only if is an edge in the graph a(t). 

Given an initial heading 6(0) = (0j(O))™ =1 and a fixed heading 6*0 in which 
agent moves at all times, for i > 0, agent z's heading evolves in discrete time 
according to the following equation: 

an 1 -n - frft) + S jeM(t)0j(*) rrn 

^ (t + 1) " TTW) ( } 

where rii(t) is the number of agents in Af{(t). 
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For a leader-following Vicsek model, we have the following correspondence of 
Theorem 13.11 

Theorem 3.2. Given a leader-following Vicsek model (V + ,a), suppose o~(oo) is 
connected. Then for any 9(0) = (6 l j(0))™ =1 and 9 , we have ]im. t ^ 00 9i(t) = 9 for 
all i = 1, 2, • • • , n. 

Proof. Note that Lemma 13.21 and Lemma 13.31 still hold for the leader-following 
case; this theorem follows from a similar argument as given for Theorem 13. II □ 

4 Conclusions and further work 

In P|, Jadbabaie et al. show that all agents shall move in the same heading if 
the neighbor-graphs are periodically jointly connected, i.e., the union of any T 
contiguous neighbor-graphs is connected for some fixed T. In the same paper, 
they also ask whether this still holds when there exists a sequence of uniformly 
bounded intervals over the discrete time such that the union of neighbor-graphs 
across each interval is connected? 

This paper has shown that all agents shall move in the same heading under a 
very weak condition that requires the neighbor-graphs to be finally jointly con- 
nected, i.e., the union of all graphs started from any time is connected. This 
result gives an affirmative answer to the question raised in [3|. 

It should be emphasized that results obtained in this paper are valid for 
many versions of the Vicsek model (or coordination multi-agent models that 
use nearest-neighbor rule to update their state) (see Remark |2.1|) . As for some 
specific versions of the Vicsek model, there have been some impressive results. 
Recently, Jadbabaie |2] has shown that, if we choose the neighborhood region to 
be open, then a necessary and sufficient condition for all headings to converge to 
the same heading is that the neighbor-graph does not change after finite steps 
and is connected. 2 Jadbabaie also notes that p. 8, last paragraph] the problem 
would be more complicated if a closed neighborhood region were chosen. It 
seems this method cannot be directly applied to other kinds of neighborhoods. 
As a matter of fact, there are often situations when agents do not have disk-like 
visibility but rather might have a cone-like field of view. 3 

Note that our results are based on the assumption that the switching signal is 
pre-specified. In our future work, we shall plan to develop a model that can ex- 
plain how the neighbor-graphs evolve over discrete time, and determine sufficient 
conditions for coordination of multi-agents in terms of these agents' initial states. 
Another thing that should be stressed is that coordination results obtained in 
ours and in [3| are all asymptotic. It will be interesting to devise other local 
updating rules such that, using these rules, coordination will be reached quickly 
and without centralized control. This work is currently being undertaken. 

2 A similar result was obtained by Savkin |B} for a simple version of Vicsek model where the 
headings 0j(i)'s take their values from a certain finite set. 

3 Such a situation has been investigated (though for continuous time) by Lin et al. 



10 



References 

[1] G. Flierl, D. Grunbaum, S. Levin, and D. Olson. From individuals to aggre- 
gations: The interplay between behavior and physics. Journal of Theoretical 
Biology, 196(4) :397-454, 1999. 

[2] A. Jadbabaie. On distributed coordination of mobile agents with changing 
nearest neighbors. Technical Report, University of Pennsylvania, Philade- 
phia, PA, 2003. 

[3] A. Jadbabaie, J. Lin, and A.S. Morse. Coordination of groups of mobile 
autonomous agents using nearest neighbor rules. IEEE Transactions on 
Automatic Control, 48(6):988-1001, 2003. 

[4] Z. Lin, M. Broucke, and B. Francis. Local control strategies for groups 
of mobile autonomous agents. IEEE Transactions on Automatic Control, 
49(4):622-628, 2004. 

[5] C.W. Reynolds. Flocks, birds, and schools: A distributed behavioral model. 
Computer Graphics, 21:25-34, 1987. 

[6] A. Savkin. Coordinated collective motion of groups of autonomous mobile 
robots: Analysis of vicsek's model. IEEE Transactions on Automatic Con- 
trol, 49(6):981-983, 2004. 

[7] J. Toner and Y.H. Tu. Long range order in a two dimensional xy model: 
How birds fly together? Physical Review Letters, 75:4326-4329, 1995. 

[8] J. Toner and Y.H. Tu. Flocks, herds, and schools: A quantitative theory of 
flocking. Physical Review E, 58(4):4828-4858, 1998. 

[9] T. Vicsek, A. Czirok, E. Ben Jacob, I. Cohen, and O. Schochet. Novel type 
of phaze transitions in a system of self-driven particles. Physical Review 
Letters, 75:1226-1229, 1995. 

[10] K. Warburton and J. Lazarus. Tendency distance models of social cohesion 
in animal groups. Journal of Theoretical Biology, 150(4) :473-488, 1991. 



11 



